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ABSTRACT 



ISO [C II] 158 /iiii, [0 I] 63 fim, and H2 9 and 17 fim observations are 
presented of the reflection nebula Ced 201, which is a photon-dominated region 
illuminated by a B9.5 star with a color temperature of 10,000 K (a cool PDR). 
In combination with ground based [C I] 609 jum, CO, 13 CO, CS and HCO + 
data, the carbon budget and physical structure of the reflection nebula are 
constrained. The obtained data set is the first one to contain all important 
cooling lines of a cool PDR, and allows a comparison to be made with classical 
PDRs. To this effect one- and three-dimensional PDR models are presented 
which incorporate the physical characteristics of the source, and are aimed at 
understanding the dominant heating processes of the cloud. The contribution 
of very small grains to the photo-electric heating rate is estimated from these 
models and used to constrain the total abundance of PAHs and small grains. 
Observations of the pure rotational H 2 lines with ISO, in particular the S(3) 
line, indicate the presence of a small amount of very warm ~ 330 K molecular 
gas. This gas cannot be accommodated by the presented models. 

subject headings: ISM: clouds - ISM: Ced 201 - ISM: molecules - molecular processes - 
reflection nebulae 
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1. Introduction 

A major problem concerning the interstellar medium and circumstellar regions 
remains the identification of the dominant gas heating source. The best current model is 
photo-electric emission by dust grains following irradiation by ultraviolet photons. The 
ejected electrons carry an excess kinetic energy of a few eV which is transferred to the gas 
by elastic collisions with H and H 2 . The efficiency of this photo-electric heating for various 
densities and illuminating radiation fields has been the subject of many investigations 
(Hollenbach & Tielens 1997; Spaans et al. 1994; and references therein). Bakes & Tielens 
(1994) presented calculations for the contribution of Very Small Grains (VSGs) and large 
molecules like polycyclic hydrocarbons (PAHs) to the total heating rate. To establish a 
firm theoretical basis for interstellar heating, one therefore has to constrain the abundance 
of these putative species. 

In this work, Infrared Space Observatory (ISO) and ground-based observations are 
presented of the reflection nebula Ced 201. Since the illuminating star has an effective 
temperature of 10,000 K, Ced 201 qualifies as a cool PDR. That is, the relatively 
soft impinging radiation field decreases the photo-electric heating rate, leading to cooler 
temperatures. Furthermore, the resulting photo-electric heating rate has larger contributions 
from grain species with low ionization potentials and negative (or small positive) charges 
(Bakes & Tielens 1994). Ced 201 is thus ideally suited to investigate the role of PAHs in 
the photo-electric heating process (Spaans et al. 1994). 

Most PDR models to date are constructed for classical PDRs, i.e. PDRs illuminated 
by O or early-type B stars. The aim of this study is to use a complete set of atomic and 
molecular diagnostics to constrain the total cooling, the density, the extinction, and the 
illuminating radiation field of a cool PDR. If the physical and chemical structure of the 
source can be determined consistent with the observations through PDR models, then 
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it is possible to estimate the contribution of VSGs and PAHs to the total photo-electric 
heating rate. ISO offers a new opportunity to probe the physics and chemistry of PDRs 
by observations of the major cooling lines [C II] 158 fim and [O I] 63 fim with the Long 
Wavelength Spectrometer (LWS) and searches for the pure rotational lines of H2 with the 
Short Wavelength Spectrometer (SWS). The latter lines provide important constraints on 
the amount of warm gas in the source, which may not be explaimed by current models 
(Timmerman et al. 1996). Searches of these LWS and SWS lines in weak, extended sources 
such as discussed here have not been possible previously with airborne observatories. 

2. Observations 

Cederblad 201 is associated with BD +69°1231, a main sequence star of spectral type 
B9.5 and an effective temperature of 10,000 K. The source is located toward a = 22:12:14 
and 5 = 70:00:11 (Epoch 1950) and lies at a distance of approximately 420 pc (Casey 
1991). In this study an extended set of observations of Ced 201 is presented. The source 
is mapped around the offset position in 12 CO 2 — > 1 and 3 — > 2 line emission. Additional 
observations of line transitions of other species are performed at the maximum of the 12 CO 
2 — > 1 emission. The most important observations are obtained with the James Clerk 
Maxwell TelescopeQ (JCMT) and ISO, and are presented in the following sections after the 
NRAO results. 



4 The James Clerk Maxwell Telescope is operated by the Joint Astronomy Centre, Hilo 
(Hawaii), on behalf of the UK Particle Physics and Astronomy Research Council, the 
Netherlands Organization for Scientific Research (NWO) and the National Research Council 
of Canada. 
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2.1. NRAO Observations 

Single dish observations of the HCO + 1-0 and CS 2-1 lines were made in June 1996 
with the NRAO 12m telescope on Kitt peak fj The adopted receiver was the 3mm SIS 
dual channel mixer. For the back end, the Hybrid Spectrometer was placed in dual channel 
mode with a resolution of 47.9 KHz (0.16 km s^ 1 ). The beam width is 63" FWHM and 
the main beam efficiency 7/mb = 0.86. Pointing is accurate to 10" in azimuth and 5" in 
elevation. A weak detection of HCO + 1 — > was found, but the CS 2 — > 1 line was not 
detected (see Table 2). 

2.2. JCMT Observations 

Observations of the 12 CO 2 -> 1 and 3 -> 2, 13 CO 2 -»• 1 and 3 -> 2, C 18 2 -> 
1, HCO + 3^2, and [C I] (492 GHz/609 jttm) line transitions were performed in three 
different wavelength bands. The front-end receivers, the telescope beam sizes, and the 
efficiencies at the various frequencies are shown in Table 1. The main beam efficiencies 7/mb 
were determined from observations of planets by the JCMT staff. All observations were 
obtained in the 125 MHz configuration, corresponding to a channel width of 0.078 MHz. 
For these high resolution spectra, the Digital Autocorrelating Spectrometer (DAS) at the 
back-end of the receiver was used. All observations were made in position switching mode, 
using an offset of 30', 45' or 90' in azimuth. In this mode, the telescope integrates at the 
source position for 30 seconds and then switches to an emission-free position on the sky. 
The difference between these two signals yields the actual source contribution, assuming 
the atmospheric conditions do not vary significantly on the time scale of a minute. Low 

5 The National Radio Astronomical Observatory is operated by Associated Universities, 
Inc., under contract with the U.S. National Science Foundation. 
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order polynomial baseline fits were adopted in the reduction. 

Table 1 

Technical Details of the JCMT 
Receiver Frequency i] MB #beam (") 

RxA2 230 GHz 0.69 ± 0.03 21 ± 2 
RxB3i 345 GHz 0.58 ± 0.02 14.1 ± 0.5 
RxC2 460 GHz 0.52 ± 0.05 10.9 ± 0.5 

Table 2 shows an overview of the line transitions observed at the offset position (0",0"), 
and Figure 1 shows the 12 CO 2 — > 1 and 3^2 line emission maps. Most line profiles can 
be fitted with a Gaussian at Vlsr = —4.9 km s^" 1 , although in several spectra a smaller peak 
at Vlsr = — 5.9 km s _1 also occurs. This study will focus on the emission at Vlsr = —4.9 
km s _1 . For Ced 201, the 12 CO line profiles do not exhibit the characteristic double-peaked 
shape of self-absorption, except for a few offset positions, located near the edge of the 
mapped regions, which are not associated with the central PDR. Therefore, it can be 
concluded that the illuminating star is located between the observer and the central cloud, 
heating the gas from the outside. 
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Table 2 

Submillimeter Observations of Ced 201 



Line 


I 




AV 


Vlsr 


r.m.s. c 




K km s" 1 


K 


km s _1 


km s _1 


K 


HCO+ 1 -> 


1 


0.15 


0.85 


-5 


0.03 


CS 2 -> 1 




< 0.03 






0.03 



12 CO 2 -> 1 


21.3 


23.2 


0.86 


-4.9 


0.53 




7.6 


6.8 


1.1 


-5.9 


0.53 


12 CO 3^2 


25.3 


33.1 


0.72 


-4.9 


0.94 




7.8 


6.8 


1.1 


-5.9 


0.94 


13 CO 2 -> 1 


5.4 


6.9 


0.73 


-4.9 


0.71 




1.7 


1.5 


1.1 


-5.9 


0.71 


13 CO 3^2 


3.9 


5.5 


0.66 


-4.9 


1.0 






<0.64 




-5.9 


0.64 


C 18 2 -> 1 




<0.09 






0.09 


HCO+ 3^2 




<0.08 






0.08 


[C I] (609 pm) 




<1.3 






1.3 



b In case of no detection, the upper limit, obtained by Hanning smoothing 
once, is shown. 

c Measured per channel of 0.078 MHz. 



2.3. ISO Observations 

Observations of the fine-structure lines of neutral oxygen [O I] and ionized carbon [C 
II], as well as the pure rotational transitions of molecular hydrogen have been performed 
by ISO. The LWS (Clegg et al. 1996) AOT02 grating mode was used to obtain spectra of 
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the [O I] (63 (ttm), [0 I] (145 /zm) and [C II] (158 //m) fine structure lines, whereas the H 2 
S(l) and S(3) transitions were obtained using the SWS (de Graauw et al. 1996) AOT02 
grating mode. The [C II] and H 2 observations show unresolved line profiles, while the [O 
I] spectra only yield upper limits, as summarized in Table 3. The reduction was performed 
with the LWS version 6.0 pipeline software and the SWS reduction package (Thi 1997, 
private communication). The detection of the S(l) and S(3) pure rotation H 2 lines indicates 
warm (T > 200 K) gas. This result will be further explored through the physical models 
constructed below. 



Table 3 
ISO Observations of Ced 201 



Line 


I 




^beam 




W m" 2 






[C II] (158 /im) 


3.80 x 10" 


-15 


90" 


[0 I] (63 /im) 


< 8 x 10- 


15 


90" 


H 2 S(l) (17 /im) 


1.10 x 10~ 


-16 


14" x 27" 


H 2 S(3) (9 /im) 


5.61 x 10" 


-17 


14" x 20" 



2.4. Comparison with Other Lines of Sight 

Before discussing the construction of specific models for Ced 201, it is good to compare 
the observed [C II] 158 /im line intensity, which accounts for 50% of the total cooling rate 
in moderate density PDRs (tin < 10 4 cm -3 ), with other lines of sight. The above value 
translates into 2 x 10~ 5 erg s" 1 cm" 2 sr" 1 , which should be compared to 2.2 x 10" 3 for W3 
(Boreiko, Betz, & Zmuidzinas 1993) and 4.6 x 10" 4 for S140 (Minchin et al. 1994). Because 
the [C II] flux scales roughly linear with the impinging UV flux in this regime, care should 
be taken in comparing these numbers. It will be shown below that the strength of the 
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incident radiation in interstellar units is G = 200 for Ced 201, roughly equal to the S140 
value, whereas W3 is characterized by G ~ 5 x 10 3 . With these numbers one finds that the 
Ced 201 [C II] line is underluminous. The remaining difference between the various lines of 
sight, and the topic of this work, is the shape of illuminating radiation field. 

Since the effective temperature of the illuminating star is 10,000 K, whereas it is 
~ 30, 000 K for S140 and W3, one would expect Ced 201 to be well described by cool PDR 
models. These exhibit a much smaller [C II] /CO ratio for the higher rotational levels of 
CO, J > 3 (Spaans et al. 1994). The reasons are that 1) the soft radiation field causes 
lower temperatures and hence smaller collisional excitation rates, and 2) the abundance 
ratio C + /CO is much larger in a soft UV field since C ionization and CO dissociation occur 
in the 912-1110 A wavelength range. Indeed, the [C II] line flux is not only small in Ced 
201, which by itself may be due to abundance, temperature, and (column) density effects, 
it is also small compared to the observed CO 3-2 line with a ratio of about 100. This is 
roughly an order of magnitude smaller than the corresponding ratios for the other lines of 
sight. The fact that the high temperature [O I] 63 /zm line, at an excitation temperature of 
about 228 K above ground, is not detected at all toward Ced 201, is consistent with these 
results if the gas temperature is smaller than 150 K throughout most of the cloud (see the 
PDR models below). 

3. A Physical Model for Ced 201 

The analysis of the data and construction of a physical model for Ced 201 will proceed 
in several steps. First, the average density, temperature and abundances of observed 
species are determined from the observations through a simple excitation model. The 
chemical abundances are then reproduced through iteration of more elaborate ID and 3D 
PDR models, which incorporate the chemical and thermal balance of the medium, and 
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the geometry of the reflection nebula. The motivation for the latter is that the spatial 
variations in the temperature and chemical abundances strongly influence the observed line 
intensities. 

The intensities of emission lines are determined by the level populations and optical 
depth along the line of sight. A one-dimensional escape probability radiative transfer code 
has been applied to determine these quantities for the detected lines. The code includes 
collisional and radiative (de-)excitation processes, and computes the level populations 
and line optical depth for a given set of physical parameters: density, temperature, 
column density, line width, and the temperature of the cosmic microwave background. 
Comparison with the observed intensities then results in a first guess for the beam-averaged 
kinetic temperature, density and column density of the observed species. This radiative 
transfer code and the procedure for analyzing the data has been described in Jansen, van 
Dishoeck & Black (1994) and Jansen (1995). In Table 4 the derived physical parameters 
are summarized. The kinetic temperature and the density of the main collision partner, 
molecular hydrogen, have been determined using the line intensities of the 12 CO and 13 CO 
transitions. The optically thick 12 CO lines were used to constrain the kinetic temperature, 
i.e. this is the temperature of the (self-shielding) molecular gas, whereas the 13 CO 3-2/2-1 
line ratio was used as the main density probe (see Figure 8a of Hogerheijde et al. 1995). 
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Table 4 

Physical Parameters and Column Densities of Ced 201 



T 


(40 ±10) K 


n(H 2 ) 


(5000 ± 1000) cm" 3 


iV(CO) 


(2 ± 1) x 10 17 cm" 2 


iV( 13 CO) 


(3 ± 1) x 10 15 cm" 2 


iV(CI) 


< 2 x 10 16 cm" 2 


JV(CII) 


(4 ± 2) x 10 17 cm- 2 


iV(C 18 0) 


< 4 x 10 13 cm -2 


iV(HCO+) 


(8 ±3) x 10 11 cm- 2 


iV(CS) 


< 1 x 10 12 cm" 2 


N(01) 


< 4 x 10 18 cm -2 



It is important to realize that the C + column density is quite sensitive to the 
temperature structure of the cloud due to its excitation temperature of 92 K. For the ID 
and 3D models we shall therefore compare to the observed flux directly. Since the more 
detailed models discussed below provide a good fit to the observed line strengths, and 
indicate column densities quite similar (within the errors) to the one above, we will use 
this knowledge and the value of the C + column density to derive now the gas phase carbon 
abundance. 

The gas phase carbon budget is dominated by 12 CO, C and C + . The sum of the derived 
column densities yields for the column density of gas phase carbon iV( 12 C) = (6 ± 2) x 10 17 
cm -2 . If it is assumed that the central gas cloud of Ced 201 is spherical, then one can 
compare the gas phase carbon column density with the column density of hydrogen nuclei 
(N-h) along the line of sight. For spherical symmetry, 7V H = nnxdiameter = (7.6±2.3) x 10 21 
cm -2 , where the diameter from the 12 CO 2-1 map is 120" at 420 pc, and the error is 
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dominated by the non-circularity of the map. Thus, one finds for the gas phase carbon 
abundance A(C) = (1.2 ±0.4) x 1CT 4 . The cosmic abundance of carbon is « 4 x 1CT 4 . Note 
here that the cosmic abundances are thought to be 30% lower than the Solar abundances 
(c.f. Meyer 1997). One finds for Ced 201 that the fraction of carbon in the gas phase is 35% 
± 10%. For diffuse clouds a gas phase carbon abundance of 1.4 x 10~ 4 is derived based on 
GHRS/HST observations (Cardelli et al. 1993), consistent with the result found here. 

If one assumes a carbon abundance of 10~ 4 , that C + is the dominant carbon bearing 
species in the atomic PDR zone, and that the ambient density is above the critical density 
of ~ 3 x 10 3 cm" 3 (LTE) as motivated by the detected CO emission, then the size of 
the emitting region and the energy rate per hydrogen nucleus required to reproduce the 
observed line flux can be estimated. Clearly, various combinations of size and energy rate 
per H atom will yield identical line luminosities. Still, in LTE one cannot increase the 
emitted flux by a large amount unless the ambient temperature is below the excitation 
temperature, in which case the line would be quite weak. Therefore, taking a size of ~ 120" 
from the CO and IRAS dust maps, one cannot make the [C II] emitting region much smaller 
than this without violating the observations, i.e. above 92 K per collision one increases 
mostly the ambient temperature and not the line luminosity. From these considerations 
it follows (to first order) that the beam filling factor of the emitting gas is roughly unity 
and the ambient temperature in the [C II] emitting region is not less than the excitation 
temperature of 92 K. 

3.1. The ID PDR Model 

Previous studies of PDRs illuminated by a (interstellar) radiation field with an effective 
temperature of about 30,000 K have been performed by Tielens & Hollenbach (1985); van 
Dishoeck & Black (1988); Burton et al. (1990); Hollenbach et al. (1991); le Bourlot et al. 
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(1992); Storzer, Stutzki, & Sternberg 1996; at various levels of sophistication. Spaans et 
al. (1994) extended PDR calculations to softer (T eff = 20,000 - 6,000 K) radiation fields 
in order to study the changes in the photo-electric heating efficiency. In this work two 
different models valid for clouds illuminated by an intense but soft radiation field, the ID 
PDR model and the 3D PDR model, will be applied to, and tested for, Ced 201. 

The one-dimensional PDR model is an extension of the diffuse cloud model of van 
Dishoeck & Black (1986) and is thoroughly described in Jansen et al. (1995, also referred 
to as the Leiden PDR model). In this model, the PDR is described as an infinite plane 
parallel slab in thermal and chemical equilibrium with a radiation field incident from one 
side. The calculation of the chemical conditions in the PDR is performed in spatial steps, 
where the radiation field at a certain point is determined by calculating the attenuation 
by dust grains, and optical depth effects in the UV absorption lines of H 2 and CO which 
lead to dissociation. Moreover, shielding of CO by overlapping H 2 absorption lines (mutual 
shielding) is also taken into account (van Dishoeck & Black 1988). The chemical network 
includes 215 different species, based on 24 elements and the isotopes of H, C and O 
(Jansen et al. 1995). The thermal balance is calculated simultaneously with the chemistry. 
The heating processes include photo-electric heating by dust grains and PAHs (Bakes & 
Tielens 1994), heating by cosmic rays (( = 5.0 x 10~ 17 s~ x ), H 2 formation heating, carbon 
photo-ionization heating, and collisional de-excitation of vibrationally excited H 2 (Tielens 
& Hollenbach 1985; Sternberg & Dalgarno 1989). Cooling is provided by spontaneous decay 
of collisionally excited fine-structure levels of C + , C and O, and rotational levels of CO. 

3.2. The 3D PDR Model 

In addition to the one-dimensional PDR model, the three-dimensional model developed 
by Spaans (1996), hereafter referred to as the 3D PDR model, is applied to Ced 201. 
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Analogous to the ID PDR model, this model calculates the thermal and chemical balance 
of the cool PDR, using the stellar radiation field, the morphology, the elemental abundances 
and the (constant) density as input parameters. Nevertheless, differences between the 
models occur, due to the geometry, radiative transfer and chemical network. The radiative 
transfer in the 3D PDR model is extended from one to three dimensions, thus enabling the 
interpretation of line intensities from geometrically more complex clouds and to take into 
account the position of the illuminating star. Furthermore, the interstellar radiation field 
(ISRF, Draine 1978) is explicitly taken into account as an isotropic background radiation 
field illuminating the side of the cloud opposite to the star. 

In contrast with the escape probability method of the ID PDR model, the 3D PDR 
model uses Monte Carlo radiative transfer (Spaans & van Langevelde 1992; Spaans 1996). 
According to this method, the PDR is divided in a large number of different cells with a 
typical size no larger than the mean free path of a photon. The physical conditions differ 
in each cell, thus influencing the excitation of atomic and molecular lines. Calculation of 
the radiative transfer proceeds by determining the stimulated emission and absorption (line 
and continuum) for propagating photon packages in each cell simultaneously. Spontaneous 
emission of line photon packages occurs in random directions. The chemical network 
incorporated in the 3D PDR model is more limited, and includes only the ~ 40 most 
important observable species. Nevertheless, the network is shown to be detailed enough to 
model the C + /C/CO transition accurately (Spaans & van Dishoeck 1997). 

4. Results 

Since the central gas cloud and the illuminating star are located in the same direction, 
the abundances as a function of visual extinction, calculated by the models, can be directly 
converted into column densities by integrating over the total extent along the line of sight. 
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The aim of the constructed models is to reproduce the observed column densities, in 
particular the gas phase carbon column density iV( 12 C), the [C II] and [O I] line intensities, 
and to determine the best fit input parameters of the ID and 3D models. 

Table 5 summarizes the input parameters resulting in the best fit ID PDR model. The 
illuminating radiation field is determined by considering the central star as a black body 
with T e ff = 10,000 K, placed at a distance of R ne buia = 0.4 pc. This corresponds to a value 
of G = 200 in units of the average interstellar radiation between 13.6 and 2 eV (Spaans et 
al. 1994), which is derived from the spectral type of the star and the infrared observations 
presented by Casey (1991) in a straightforward manner. The density of hydrogen nuclei 
in Ccd 201 is found to be n H = 2n(H 2 ) + n(H) = 1.2 x 10 4 cm -3 and reproduces well 
the observed chemical abundances. The best fit column density of molecular hydrogen 
iV(H 2 ) = 1.6 x 10 21 cm -2 is a measure of the physical size of the cloud. When compared 
with the spherically symmetric case, iV(H 2 ) S ph. sy mm. = 3.8 x 10 21 cm~ 2 , its value indicates 
that the PDR is flattened along the line of sight. Note that this also implies that the 
abundance of gas phase carbon must be higher than 35%. In fact, the flattened geometry 
yields 80% for the abundance of carbon in the gas phase relative to the cosmic abundance. 
This revised gas phase abundance is thus a factor of two higher than the value of 1.4 x 10~ 4 
found for diffuse clouds. The best fit model also constrains the gas phase abundances S of 
several other elements, which are shown in Table 5. The gas phase fraction of nitrogen is 
rather arbitrarily chosen, since no observations of nitrogen bearing species were obtained. 
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Table 5 

Input Parameters for the Best 



Fit of the ID PDR Model 



-^nebula 


U.4 pc 


n H 


1.2 x 10 4 cm" 3 


iV(H 2 ) 


1.6 x 10 21 cm" 2 


5C 


0.80 


50 


0.80 


5S 


0.01 


8N 


0.70 


<5metals 


0.02 



The 3D PDR model is applied consistently with the ID PDR model. That is, the only 
input parameter of the 3D PDR model which can be modified for a constant density cloud 
is the three-dimensional geometry of the cool PDR. The best reproductions of the observed 
column densities are obtained for an oblate ellipsoid with an axis ratio of 2.4. The ID and 
3D models which yield the best fit to the observations are presented in Table 6. 
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Table 6 

Comparison Between the Observed Column Densities, 
the ID PDR Model and the 3D PDR Model 



Species 


Observed 


ID PDR 


3D PDR Model 




cm -2 


cm -2 


cm -2 


T T 

H 2 




1.6 x 10 


1.5 x 10 21 


12 CO 


(2±1) x 10 17 


2.3 x 10 17 


2.1 x 10 17 


13 CO 


(3± 1) x 10 15 


4.4 x 10 15 


3.1 x 10 15 


c 18 o 


< 4 x 10 13 


7.2 x 10 13 


3.1 x 10 13 


c 


< 2 x 10 16 


5.2 x 10 16 


3.2 x 10 16 


c+ 


(4 ±2) x 10 17 


3.2 x 10 17 


3.7 x 10 17 


cs 


< 1 x 10 12 


6.3 x 10 11 


5.1 x 10 11 


HCO+ 


(8 ±3) x 10 11 


4.1 x 10 11 


6.8 x 10 11 



4.1. Chemical Structure 

The ID PDR model yields the abundances as functions of depth for approximately 
200 species, of which the most important ones are plotted in Figure 2. The assumed 
isotope ratios are [ 12 C]/[ 13 C] = 60 and [ 16 0]/[ 18 0] = 500. At the illuminated side of Ced 
201, ionized carbon is the dominant carbon bearing species. Due to line shielding and 
attenuation by intervening dust, the ionized carbon is converted into neutral carbon and 
carbon monoxide at Ay ~ 1 mag (Tielens & Hollenbach, 1985, Black & van Dishoeck, 1987, 
and references therein). Nowhere in the cloud does neutral carbon become the dominant 
carbon bearing species. 

The abundances of several species determined by the 3D PDR model are shown in 
Figure 3. Consistent with the ID PDR model, the C + /C/CO transition occurs at an 
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extinction of Ay ~ 1 mag. Due to the ISRF which illuminates the cloud from all directions, 
a second C+/C/CO transition zone is present at Ay ~ 2 mag on the far side of the PDR. 
Atomic hydrogen is the dominant form of hydrogen at the illuminated edge, but the 
shielding of H 2 is very efficient and leads to a sharp transition. 

4.2. Thermal Balance and Small Grains 

Both the 3D PDR model and the ID PDR model compute the thermal balance in Ced 
201. Figures 4 and 5 show the heating and cooling rates as functions of depth determined 
by the ID PDR model. Except for the heating by cosmic rays, all heating mechanisms 
decrease with depth into the cloud. Photo-electric emissions by PAHs and dust grains are 
the dominant heating mechanisms in Ced 201. The heating rates of these two processes are 
equal, although only 10 % of the solid state carbon is incorporated into PAHs in the best fit 
model. This shows how quantum effects in PAHs can contribute significantly to the heating 
process in cool radiation fields (Bakes & Tielens, 1994). Cooling is provided by radiative 
decay of collisionally excited species. The cooling rates due to the main coolants are shown 
in Figure 5, and depend strongly on the abundances of these species. The observations do 
not provide good constraints on the column density of atomic oxygen, and a gas phase 
fraction of 3.2 x 10~ 4 is therefore assumed (Meyer 1997). 

The thermal balance is established at an equilibrium temperature profile determined 
by the heating and cooling rates. Figures 6 and 7 show these profiles for the ID PDR 
model and the 3D PDR model, respectively. The temperature profiles calculated by both 
models are approximately the same. The temperature at the illuminated edge is 168 K 
according to the 3D PDR model, and 169 K for the ID PDR model. The temperature 
decreases with increasing depth, and is for both models ~ 40 K at Ay = 1 mag, where the 
temperature tracer carbon monoxide becomes abundant. According to the ID PDR model, 
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the temperature gradually declines at the dark side of the PDR. Due to the ISRF however, 
the 3D PDR model predicts that the temperature will increase again, after reaching its 
minimum of 22 K at Ay = 1.8 mag, and becomes 34 K at the dark side of Ced 201. 

For the best fit models to the chemical abundances, one can determine the model 
dependent, i.e. with the distance to the star fixed by the infrared observations, contribution 
to the heating rate from PAHs and VSGs. The models indicate that 10% of the gas phase 
carbon is locked up in the form of PAHs, i.e. -Dpah = 0.1, where PAHs are defined as any 
carbonaceous particle of linear size less than 15A. Additional models were run without any 
contribution of PAHs and VSGs to check the robustness of this result. The resulting thermal 
balance in the latter case yields temperatures at the edge of the cloud which are lower by 
50%, and strongly underestimates the observed [C II] flux. The fraction -Dpah = 0.1 also 
yields good agreement with the observed CO excitation temperatures. The precise value 
for -Dpah depends on the model. Nevertheless, a firm conclusion is that heating by PAHs 
and VSGs is required to explain the observed properties of Ced 201, provided no heating 
sources other than photo-electric heating and cosmic rays are important. 

A final crucial comparison comes from the predicted [C II] 158 [im and [O I] 63 [im 
line intensities, because they reflect the thermal balance of the medium most strongly. In 
the ID case one finds values of 1.1 x 10~ 5 for [C II] and 1.9 x 10 -5 erg s _1 cm~ 2 sr _1 for 
[O I]. The 3D results are systematically higher and yield 1.8 x 10~ 5 and 2.7 x 10~ 5 erg s _1 
cm~ 2 sr -1 , respectively. Comparison with the observed [CII] intensity would thus favor the 
3D model, but the difference is small. 
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4.3. Pure Rotational H 2 Lines 

The H 2 J=3 and J=5 levels in the vibrational ground state lie at ~ 1,000 K and 
~ 2, 500 K above ground, and therefore trace warm molecular gas. The ISO observations of 
the H 2 S(l) and S(3) pure rotational transitions are not included in the PDR models, since 
at kinetic temperatures of T < 100 K, i.e. in a large part of Ced 201, the absolute intensities 
of these lines are negligible. For higher temperatures these lines are useful temperature 
tracers because they are optically thin and the level populations are in thermal equilibrium. 
The ratio between the observed transitions is 7(9/xm)//(17/xm) = 0.69, which indicates a 
kinetic temperature of T ~ 330 K (Black & van Dishoeck 1987). 

This observational result is clearly at variance with the model predictions. Particularly 
since the effective temperature of the illuminating is low. It appears that an additional 
heating source is required in some part of the cloud. If we assume that a thin layer of 
hot gas is present at the illuminated edge or any where else inside the cool PDR, then it 
follows from the absolute intensities that this layer is not thicker than 5A y ~ 0.05 mag for 
a mean H 2 fraction of 0.2. We therefore feel that the value of -Dpah remains secure since 
the chemical constraints which enter into it, derive from regions in the cloud with Ay > 0.3 
mag. The upper limit of 5 x 10~ 5 erg s _1 cm -1 sr" 1 for the [O I] 63 /im line and its density 
dependence of 1.0 x 10~ 8 nH 2 erg s -1 cm -1 sr _1 for temperatures higher than 228 K and 
subthermal excitation, indicate that an upper limit for the H 2 gas density of 5 x 10 3 cm~ 3 
when a unity beam filling factor is adopted. /,From the CO maps we estimate a beam 
dilution factor of not more than two, so the upper limit is consistent with the CO density 
estimate and vibrational de-excitation heating by H 2 can be ruled out. 

Furthermore, the continuum intensity observed by IRAS is around 2 x 10~ 3 erg s _1 
cm -2 sr -1 , whereas the [C II] intensity is 2 x 10~ 5 erg s" 1 cm" 2 sr -1 . This translates into a 
heating efficiency for the gas of 1%, a generic value for PDRs. The temperature of 300 K 
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itself requires a G of 5 x 10 3 for an ambient H 2 density of 5 x 10 3 cm -3 , incompatable with 
the value Go — 200 from the illuminating star. All in all, the high temperature molecular 
gas remains to be explained, and does not seem to be accommodated by the process 
of photo-electric heating. In fact, observations presented by Witt et al. (1987) on the 
scattering properties of dust in Ced 201 indicate that there is a narrow size distribution of 
grains skewed toward larger than average particle sizes. This renders a larger contribution 
to the photo-electric heating by small grains unlikely. 

One might argue that the H 2 levels are populated through UV pumping. Even 
though the reflection nebula is relatively close to the illuminating star, the low effective 
stellar temperature strongly quenches the 912-1110 A flux relative for H 2 fluorescence. A 
straightforward calculation shows that UV irradiation can account for at most 20% of the 
S(3) line, thereby rendering it a minor contribution. 

An additional heating source like turbulent dissipation could be present in the cool 
PDR (Falgarone & Puget 1995). This requires the input of kinetic energy on the scale of 
the cloud, possibly through a weak (3-5 km s -1 ) shock. A C-shock of 7 km s -1 into a gas 
of pre-shock density 10 4 cm" 3 would also suffice, but would likely lead to an observable 
[O I] 63 /iiii line. Finally, the case of Ced 201 appears not to be unique. Observations 
obtained by Timmermann et al. (1996) for S140 indicate similarly warm gas. It would be 
quite interesting if more of such warm regions show up in H 2 rotational line data. 



5. Conclusions 

The ID PDR and the 3D PDR model reproduce the chemical abundances and the 
thermal balance of the cool PDR Ced 201 as derived from the observations. All observed 
lines, obtained with the JCMT, ISO and NRAO, are accommodated by these models, except 
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for the pure rotational transitions of H 2 . The observed H 2 S(l) and S(3) line intensities 
indicate the presence of hot gas with a kinetic temperature of T ~ 330 K, which is probably 
located in a thin layer of 5 Ay ~ 0.05 mag at the illuminated edge of the cool PDR. The 
best fit models indicate that the gas phase carbon abundance is 50% of Solar, and that 10% 
of the available carbon atoms is in the form of PAHs and VSGs. 
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Fig. 1. — Contourplot of Ced 201. In the left panel the integrated line intensities with 
spacings of 60" of the 12 CO 2 — >• 1 transition are shown. The contours represent values of 
0.6, 5, 10, 15, 20 and 25 K km s _1 . The right panel presents the integrated line intensities 
of the 13 CO 2 — >• 1 line with spacings of 20". The contour values correspond to 2, 4 and 6 K 
km s -1 . The 13 CO 2 — »• 1 emission is also indicated in white in the left panel. 

Fig. 2. — Abundances of carbon bearing species in Ced 201 as functions of visual extinction 
into the cloud, for the ID PDR model. 

Fig. 3. — Abundances of several species as functions of visual extinction, for the 3D model. 

Fig. 4. — The most important heating rates as functions of visual extinction (ID PDR 
model). 

Fig. 5. — Cooling rates originating from the main coolants as functions of visual extinction 
(ID PDR model). 

Fig. 6. — Temperature distribution throughout the cloud, according to the ID PDR model. 

Fig. 7. — The kinetic gas temperature as a function of extinction, according to the 3D PDR 
model. 
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